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Abstract. We explain the exact solution of the 1+1 dimensional Kardar-Parisi- 
Zhang equation with sharp wedge initial conditions. Thereby it is confirmed that 
the continuum model belongs to the KPZ universality class, not only as regards 
to scaling exponents but also as regards to the full probability distribution of the 
height in the long time limit. 



1 Introduction 



Over the past year there has been significant advances in the understanding of 
surface growth in 1+1 dimensions, both on the experimental and theoretical side. 
The experiments on droplet growth for a thin film of turbulent liquid crystal were 
presented at StatPhys24 by K.A. Takeuchi [TJ. We will report here on the theoretical 
findings, but use the occasion to put the results in a wider context. We also take 
the liberty to comment on the experimental situation. 

In their seminal 1986 paper [2] Kardar, Parisi, and Zhang (KPZ) proposed an 
evolution equation for growing interfaces which in 1+1 dimensions, of interest here, 
reads 

|*=iA(|*)' + ^ + VS,. (1) 

(0Q) is a stochastic evolution equation for the height function h(x, t), hence has to be 
supplemented by the appropriate initial conditions. The height h depends on the 
spatial location x G R and on time ( e K + . ^X(dh/dx) 2 is the nonlinear growth 
velocity depending quadratically on the local slope. The Laplacian smoothens the 
height profile with relaxation coefficient v > 0. r](x,t) is space-time white noise, i.e. 
T) is Gaussian with mean zero and covariance {rj(x,t)r)(x',t')) = 5(x — x')5(t — t'). 
\f~D is the intensity of the noise and y/Drj models the random nucleation events at 
the interface. We consider only the infinite line without boundary conditions. Finite 
geometry and possibly boundary sources are of physical interest. The corresponding 
large deviations have been studied by Brunet and Derrida [3]. In this article our 
main focus is on typical fluctuations. Their finite volume corrections have been little 
explored so far. 

Already in their original contribution, KPZ noted that (CQ) can be transformed to 
a linear equation at the expense of turning the additive noise rj into multiplicative 
noise. More precisely, one introduces the Cole-Hopf transform 

Z(x, t) = exp [(X/2u)h{x, t)] . (2) 

Then Z satisfies 

^Z{x, t) = v-^Z(x, t) + (\v / D/2v)r l (x, t)Z{x, t) . (3) 

To "solve" ([3]) we introduce the auxiliary standard Brownian motion b(t), t > 0, 
with variance t. Then, by the Feynman-Kac formula, 

Z(x,t) = E^exp [a dsr)(b(s), s)] exp [(\/2v)h(b(2vt), 0)] ) (4) 

with a = (2v)~ 3 / 2 \D 1 / 2 . To distinguish from the white noise average denoted by 
(■), the average over the Brownian motion starting at x is denoted here by E x (-). 
h(x, 0) is the initial height profile. According to ([2]) the polymer has inverse stiffness 
z/, which for notational convenience has been absorbed in the time scale of b(t). 
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Z(x,t) can be viewed as the random partition function of the directed polymer 
b{s), < s < 2ut. The polymer is subject to the random potential r\ with strength 
a. To compute the potential energy of the polymer one simply has to sum rj along 
the location of the polymer chain. Inverting (j2J), the height is the random free energy 
of the polymer. Thus the KPZ equation is equivalent to a model in the equilibrium 
statistical mechanics of disordered systems. Since the free energy valleys are in 
two-dimensional space-time, the directed polymer is easier to visualize than the 
high-dimensional configuration space of Ising spin glasses. 

We remark that the KPZ equation and the Cole-Hopf transformation generalize 
in the obvious way from 1 + 1 to d + 1 dimensions. 

There are books, comprehensive reviews, and a large number of original articles 
on growth models in the KPZ class jH El U\- Of course, we cannot provide here a 
comprehensive summary. But very schematically, it seems reasonable to distinguish 
two in essence disjoint lines of research. 

• Scaling theory, critical exponents, and phase diagram in general dimension, start- 
ing with KPZ in 1986, 

• exact probability density functions (pdf) for growth models in 1 + 1 dimensions, 
starting with Johansson [8] in 2000. 

We comment on a few general features. 

Scaling theory, phase diagram. Interface motion in the KPZ class is characterized by 
local growth rules with single growth events being statistically independent in space- 
time. There is no mass transport along the interface. The respective bordering bulk 
phases are statistically homogeneous and non-critical. As a rule, in the course of 
time the interface develops a statistically self-similar structure, when viewed on suf- 
ficiently coarse scales. The self-similar structure has been amply confirmed, mostly 
by Monte-Carlo simulations of various discretized models, for which there are many 
options: space can be discretized to a lattice Z, the height may take only integer 
values, and one can implement updates in discrete time steps. Slightly more distant 
from (CQ) would be models as the off-lattice Eden growth model, where hard balls of 
equal size are added at random one by one to the current cluster of balls [9]. In fact, 
convincing numerical solutions of the 1 + 1 dimensional KPZ equation are only fairly 
recent [10]. All these models constitute the KPZ universality class. One expects 
that in the infrared scaling, i.e. for large distances and long times, essentially all 
microscopic details become irrelevant and universal statistical laws emerge. 

Amongst the many advances the most important findings can be summarized as 
follows. 

(i) Phase diagram. For dimension d < 2 there is only a strong coupling phase, ex- 
cluding exceptional points where the effective coupling constant vanishes. For d > 2 
there is a weak coupling and strong coupling phase. In the weak coupling phase 
the nonlinearity is irrelevant and fluctuations are according to the linear (A = 0) 
Edwards- Wilkinson equation. Over the last ten years probabilists have proved in 
detail many of these claims. They mostly consider a directed polymer, where b(t) is 
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replaced by a discrete time random walk and r)(x,t) by a collection of independent 
random potentials V(j,t), j 6Z, t £ Z + . In particular, it is proved that in the weak 
coupling regime the directed polymer has the statistics of a Brownian motion with 
a diffusion constant computed from the annealed version, see [IT] for an overview 
and [12] for recent progress. In fact, the bounds on the phase diagram are fairly 
sharp. From the perspective of spin glasses a natural quantity to consider is the 
ratio (Z(x, t) 2 )/(Z(x, t)) 2 . For d > 2, at small coupling the ratio stays bounded, 

(2) 

while beyond a critical coupling, Ac , it diverges exponentially. Based on theoreti- 
cal arguments and numerical simulations it has been claimed that the true critical 
coupling A c = Ac [13] - Rigorous bounds establish A c > A c , although with a fairly 
small difference [T4"| [15] . 

(ii) Critical exponents. There is a static critical exponent, \i which governs the 
height-height correlations in the steady state at fixed slope and a dynamical critical 
exponent, z. Both are related through 

z + X = 2. (5) 

For d = 1, in the steady state of the KPZ equation the slope dh/dx is distributed 
according to spatial white noise, hence X = \ and z = 3/2. For d > 1, extensive 
numerical simulations are available. In the analysis of the numerical data one has 
to assume, implicitly or explicitly, a particular form of the finite size scaling, for 
which no systematic theory is available. This is one essential limiting factor in the 
attempt to compare with theories making definite predictions for the value of the 
exponents. 

(Hi) Nonuniversal constants. All microscopic details of either an experimental re- 
alization or of a particular theoretical model are subsumed in three nonuniversal 
constants [T6J H7j . One considers a particular direction along which the interface 
has a non-random time-independent macroscopic slope u. Then the first coefficient 
is simply the steady state growth velocity reached for long times. The second one 
is the effective coupling constant, defined by the Hessian (<9 2 / ' duidu^v^. The KPZ 
theory assumes that the eigenvalues of the Hessian are either all strictly positive or 
all strictly negative. (If the eigenvalues have a different signature the fluctuations 
will be Gaussian [13 [19].) The third coefficient, A, is the amplitude of the steady 
state correlation function at slope u, i.e. ({h(x) — h(x')) 2 ) = A\x — x'\ 2x for large 
| a; — x'\. In general, all nonuniversal coefficients depend on the slope u. 

In the present context, the central unsolved problem is the upper critical dimen- 
sion. One prediction, based on asymptotic expansions, is that z increases with d 
and reaches the asymptotic value 2 as d — > oo. Other field theory based approaches 
yield the upper critical dimension d — 4. For d > 4 one finds z — 2, like in the 
weak coupling phase, but there are other statistical properties which still distinguish 
between weak and strong coupling phase. We refer to [201 EU [22] and the recent 
letter by Canet et al. [23] with references to earlier literature. <0> 

Exact solutions in 1+1 dimensions. Ulam proposed the following combinatorial 
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problem: One considers a random permutation p(l) . . .p(N) of 1 ... N. Any given 
permutation will have increasing subsequences, e.g. 641 523 has the increasing 
subsequences 4 5, 1 2 3 and others. Ulam asked for the random length, £ N , of the 
longest increasing subsequence (there could be several ones with the same length). 
One finds that typically £n = 0{s/N) for large N. So the more detailed question 
deals with the fluctuations relative to y/~N. In a famous contribution, Baik, Deift, 
and Johansson [24j established that the fluctuations are of order iV 1//6 . While the 
exponent has been anticipated from numerical simulations, the real surprise was 
that the pdf of the fluctuations turned out to be identical to the Tracy- Widom dis- 
tribution known from the Gaussian unitary ensemble (GUE), thus establishing a 
completely unexpected link to random matrix theory. Johansson [8] extended these 
results to the single step growth model (alias TASEP, the totally asymmetric simple 
exclusion process). In this model space is discrete, j G Z, and a height variable h(j, t) 
takes only integer values satisfying the single step constraint \h(j+l,t) — h(j,t)\ = 1. 
In the stochastic update, independently at each local minimum the height is in- 
creased by 2 with probability dt and stays put with probability 1 — dt. As initial 
condition we set h(j, 0) = Johansson proved that for the height at origin it holds 

h(0,t) = ±t + 2- 1 /V/ 3 £tw (6) 

for large t. Here £tw is a GUE Tracy- Widom distributed random variable, as will 
be discussed further below. As predicted by the KPZ theory the fluctuations of the 
height are of order t 1//3 . Also, the numerical coefficients in (p]) are the nonuniversal 
constants explained in (Hi) above. 

In [25] it was shown that Ulam's problem is equivalent to the polynuclear growth 
model (PNG) with droplet geometry, i.e. the height function has typically the shape 
of a semicircle. As a surprising consequence, it was realized that while the size of the 
fluctuations is always of order i 1 / 3 , the statistics still carries some information on 
the initial data [251 [261 W\ I2HI W\ ■ E.g., in case of the PNG model with flat initial 
conditions one finds that ([6]) is still valid, but the pdf of the random amplitude is 
given by the Tracy- Widom distribution of the Gaussian orthogonal ensemble (GOE). 

The discoveries from 2000 initiated ongoing activities. We refer to a book and 
reviews [30[ EH [321 [33], where more details can be found. <0 

With the exact solutions obtained since 2000, one has now available a more strin- 
gent test for a one-dimensional growth model to be in the KPZ universality class: 
Not only must be the size of the fluctuations of order t 1//3 , but, for example in the 
droplet geometry, the random amplitude must be GUE Tracy- Widom distributed. 
For short, let us call this property the universal one-point pdf. One would expect 
that any model in the KPZ class shares the universal one-point pdf. In fact, one 
should regard it as part of the definition for a model to be in the KPZ universal- 
ity class. So far the evidence comes from the few models which allow for exact 
computations. One might have guessed some activity from the side of Monte-Carlo 
simulations. In the early period the statistical sampling was limited. But cumulants 
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up to the fourth order [T6l [T7] and even the full height distribution function [34] 
have been computed for various models, in retrospect achieving surprisingly good 
agreement with the exact solution in 1+1 dimensions. More recently large scale 
simulations have been carried out, e.g. for off-lattice Eden models in the plane. 
Exponents are measured with high precision [9], but the full pdf has not yet been 
determined. 

The list of exactly soluble models is short. 

• PNG model, TASEP with discrete and continuous time update, PASEP. These 
models are unified in the directed polymer picture. They are zero temperature ver- 
sions, i. e. one studies the fluctuations of the ground state energy of the directed 
polymer. In the discrete time TASEP the directed polymer lives on the lattice Z 2 
with the time direction along (1,1). The droplet geometry means that both end- 
points of the polymer chain are fixed. The random potential V(i,j), (i,j) G Z 2 , 
is independent, identically distributed with a one-sided geometric distribution, i.e. 
Piob({V(i, j) = n}) = (1 — a)a' n ' for n < and otherwise with < a < 1. The 
PNG model corresponds to the Poisson limit of rare events, a < 1, while the con- 
tinuous time TASEP corresponds to the limit of a one-sided exponential. 

• The partially asymmetric simple exclusion process (PASEP). This is the single step 
model with a modified rule for updates. As before local minima are increased by 2 
at rate p and, in addition, local maxima are deceased by 2 with rate q, q + p = 1. 
The symmetric case p = q = \ has Gaussian fluctuations. For p > q, Tracy and 
Widom [35J prove the asymptotics of the form ([6]). 

• KPZ equation. The KPZ equation (DO) has the static exponent x = \i as a conse- 
quence of the exactly known steady state, and hence the dynamic exponent z = 3/2. 
To establish the universal one-point pdf requires a detailed computation, as will be 
explained in the remainder of this contribution. For convenience we set A > 0. The 
initial conditions are the sharp wedge h(x, 0) = —\x\/8 in the limit 5 — > 0. Then, 
for large t, 

h{0,t) = -^( 7t ) 3 + 21oga + 7^TW, (7) 

where j t = (aVt) 1 / 3 . ([7]) should be compared with ([6]). The prefactors are the 
nonuniversal constants specific for the KPZ equation, while the logarithmic shift 
has its origin in the slightly singular initial condition. In fact, ((7]) remains valid 
for every t provided (tw is replaced by the random variable C,f Its variance stays 
bounded as a function of t and its precise pdf is the main result of our contribution. 

To close the introduction, we outline some experimental activities. 

(i) Ballistic deposition |36j. A rough surface is generated by raining material onto 
a solid substrate. The randomness results from the incident beam. One difficulty is 
the suppression of surface diffusion. The incident velocity must be small to ensure 
a proper local attachment. Also the intensity must be small so to deposit on a re- 
laxed surface. The KPZ fluctuations become visible only in depositing many layers. 
The surface roughness is measured through small angle scattering of X-rays. With 
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these severe constraints, at best a qualitative confirmation of the KPZ roughness is 
obtained. 

(ii) Smoldering paper [37] . One fixes a piece of paper in a frame and, after suitable 
chemical preparation, smolders it from the bottom. It is argued that the one- 
dimensional combustion front is governed by the KPZ equation. The random noise 
originates from the intrinsic structure of the paper. Clearly, the geometry corre- 
sponds to flat initial conditions and one would like to observe the 1/3 exponent and 
the GOE Tracy- Widom distribution. To discern such fine statistical properties the 
experiment has to be repeated many times with identical KPZ parameters, which 
is difficult to achieve. A total of 18, resp. 21, successful burns are recorded for 
horizontal, resp. tilted ignition. Using two fitting parameters for the probability 
density function, reasonable agreement with the theory is obtained. 

(iii) Facet edge fluctuations [3H]. The experiment is based on the observation that 
the 1+1 dimensional KPZ equation governs also the thermal fluctuations of the edge 
of a crystal facet (which is a one-dimensional structure). To measure directly the 
spatial edge fluctuations is essentially impossible. However one can observe the dy- 
namical fluctuations of the facet edge which via detailed balance are related to the 
static fluctuation exponent. The experiment clearly rules out random walk like fluc- 
tuations of the facet edge, characteristic for an isolated step, and favors the reduced 
KPZ roughness. Qualitatively the reduction is caused by the constraint through the 
neighboring step away from the facet and by entropic constraints for deviations into 
the facet. 

(iv) Turbulent liquid crystal pQ. One prepares a thin film of liquid crystal with 
dimensions 16mm x 16mm x 12/mi. The film is confined by two transparent elec- 
trodes and subject to oscillating electric fields. A particular point, temperature 
25 °C and voltage 26 V at 250 Hz, in the nonequilibrium phase diagram is chosen 
with care, such that the so-called DSM 2 phase is stable, while the DSM 1 phase is 
unstable. However the time for spontaneously nucleating the DSM 2 phase in the 
bulk DSM 1 phase is long compared with the time of a single experimental run. 
Also unavoidable nucleations at the border of the cell are fairly unfrequent. The 
system is prepared in the homogeneous DSM 1 phase. A laser pulse plants a seed 
of the DSM 2 phase which then grows within 35 sec to a droplet filling the region 
seen by the camera. This is precisely the geometry used in the exact solution of the 
KPZ equation. To have flat initial conditions the laser pulse plants a line seed. The 
experiment is repeated of the order of 1000 times. The DSM 1 phase is turbulent, 
hence good mixing in the region away from the droplet is ensured. For a single probe 
the physical parameters remain constant over long periods. By spatial isotropy the 
average shape of the droplet is a disk which improves substantially the analysis of 
the statistical data. First the nonuniversal constants are determined with good pre- 
cision. Thereby the scales are fixed and no fitting parameters are used. The scaling 
exponents, the GUE and GOE Tracy- Widom pdf, and the spatial height-height cor- 
relations are measured. Very good agreement with theory is achieved. Also finite 
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time corrections are determined. 



2 The Cole-Hopf solution and its approximations 

As written, the KPZ equation is not well-defined mathematically. The noise is 
singular and even if one takes the smoothening by the Laplacian into account the 
solutions to the linear part of the equation is rough. Thus the nonlinearity requires to 
multiply pointwise two distributions which is an ambiguous mathematical operation. 
The proper definition of solutions to the KPZ equation is ongoing research. For our 
purposes the Cole-Hopf transform suffices. Of course, it still reflects the difficulties 
mentioned: the action integral is not well-defined for typical realizations of Brownian 
motion and of white noise. Fortunately, as in other two-dimensional field theories, an 
infinite energy renormalization suffices to make sense out of (j3J). Thus the strategy is 
to suitably approximate Z(x,t) and then to define the solution to the KPZ equation 
as 

h{x,t) = {2u /X) log Z(x,t) . (8) 

We will study exclusively the sharp wedge initial conditions, which under Cole- 
Hopf translates to 

Z(x,0) = 6(x). (9) 

Currently there are five different approximation schemes, which all yield the 
same limit stochastic process. We discuss each one separately. 

(A) Multiple Ito integrals. We expand the exponential of fll]) and average over 
Brownian motion which yields as n-th coefficient 

a n dti...dt n dx 1 ...dx n T\r}(x j7 t j )p (x 1 ,ti,...,x n ,t n ,x,2vt) 

Ja<t 1 ...<t n <2ut JR n =1 

(10) 

with po the joint probability density for the Brownian motion to start at 0, to be 
at Xj at time tj, j = 1, . . . , n, and to be at x at time 2z/t. We now integrate in x 
and interpret the remaining time integrations as multiple Ito integral in the forward 
discretization. Then by Ito's lemma 
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(Z(x,t) 2 } = J2 a2n / dtx...dt r , 

n=0 J0<t 1 ...<t n <2ut 



x / dx! . . .dx n po(x!,ti, . . . ,x n ,t n ,x,2vt) 2 < oo , (11) 

since the n-fold integral can be bounded by C/fal) 1 / 2 . Hence Z(x,t) is a well 
defined random variable. In fact, it can be shown that Z(x,t) > and that Z(x,t) 
is continuous in x and in t > with probability one [39J. 

(B) Colored noise. One introduces the mollifier (p K (x) = K,ip(nx), with (p > 0, (p of 
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rapid decrease, (f(x) = tp(—x), J dx(p(x) = 1, and spatially smears the white noise 
to T] K (x,t) = J dx'ip K (x — x')r](x',t). The Cole-Hopf transformation remains valid 
and the action is now a properly defined integral. The white noise average over the 
partition function is given by 
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(Z K {x,t)) = (E (exp a / dsrj K (b(s), s) 5{b(2ut) - x)) ) 

L Jo J 

= p (x, 2vt) exp[i«V K * <Pn(0)t] , (12) 



where * denotes convolution and po(x,t) = (2nt)~ 1 ^ 2 exp[— x 2 /2t] is the transition 
probability of standard Brownian motion. As k — > oo one has r\ K — > r), but (Z K (x, £)) 
diverges as exp[^a 2 ip*ip(0)Kt]. On the level of the KPZ equation this means that the 
average velocity of the interface, v K = \a 2 ip K * <f K (0), diverges linearly in k. But it 
can be proved that in the frame moving with constant velocity v K in the /i-direction 
one nevertheless has a well-defined limit as k — > oo. In other words 

lim p (x, t)Z K (x, t)/ (Z K (x, t)) = Z(x, t) (13) 
with Z(x,t) as defined in item (A). 

(C) Lattice directed polymer. The integral (j3J) defining the random partition func- 
tion is discretized. One replaces the Brownian motion b(t) by the random walk u 
on the lattice Z 2 . uj starts at (0, 0) and moves with probability | either up or right. 
At each site of the lattice there is independently a unit Gaussian random potential 
f)(i,j). Then the discrete approximation to (j4j) reads 

Z N (x,t)= J2 2- tN e~^\ (14) 

ui:(0,0)~>{(t+x)N,(t-x)N) 

where the sum is only over paths with endpoint ((t + x)N, (t — x)N), \x\ < t, and 
the energy of the walk u is defined through 

E{u)= J2 (15) 
(i,i)eoj 

As proved in [3D], in the limit — > oo, (3 4 N = a fixed, 

lim Z N (x,t) = Z(x,t) (16) 

iV->oo 

as a stochastic process and Z(x, t) as in item (A) with parameters v = | and D = 1. 

(D) The attractive 5-Bose gas. For integer moments of the partition function with 
colored noise as in (B) one can carry out the Gaussian average over r\ with the result 

{Z K (x,t) n } = (0|e- 2 ^ (K) |x). (17) 
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On the right one has a matrix element of a quantum propagator for n particles on 
R. More explicitly 

n 

\a 2 ^2<p K *<p K (xi-Xj) (18) 

with Xj G IR the position of the j-th quantum particle and \x) denoting the quantum 
state where all n particles k — y oo, Z K {x,t)/(Z K {x,t)) converges to 

Z(x,t)/p (x,t) with Z(x,t) as in (A), while on the right one finds the hamiltonian 
for n quantum particles with an attractive 5-potential and the self-term omitted. 
Therefore 

(Z(x,t) n } = (0\e- 2utH "\x), (19) 

n rj2 n 

ff --Eb-KE^- s i)' ( 20 ) 

3=1 ' <#i=l 

While (TT9"j) . together with (I2"0~j) . is a correct identity, it cannot be used to define the 
pdf of Z(x, t), since log(Z(x, t) n ) = n 3 and the moment problem defined by (fT9]) has 
many solutions [41] . 

(E) Weak ASEP. We consider the single step model with asymmetry q, q > |. The 
PASEP occupation variables are denoted by rjj(t), j G Z, rjj(t) = 0, 1. The initial 
height profile is h(j, 0) = For this initial condition the height at time t is given 
by 

j 

h(j,t) = -2 Vt(t)+j. (21) 

e=-oo 

We have seen that for the directed polymer the KPZ equation corresponds to a 
weak noise approximation. In a similar spirit, the PASEP should approximate the 
KPZ for weak asymmetry. The asymmetry must be carefully chosen. If it is too 
weak, one arrives at the Gaussian theory with A = and, if the asymmetry is 
too strong one misses the approximation through a stochastic partial differential 
equation altogether. The correct choice is as follows: first we scale space-time 
diffusively, i.e. j = L^ 1 ^] with |_-J denoting integer part and time as e~ 2 t, t = 0(1). 
< e 1 is the scaling parameter. With this choice we set q + p = 1 and 
q — p = p>/e, (3 > 0. To recall the e-dependence of q, we write h £ (j,t) instead of 
h(j,t) in the definition (l2"Tj) . As established in [121125], as convergence of stochastic 
processes it holds 

lim ^Ph £ { [e~ x x\ , e-H) + i^te" 1 - ±(3H - log(2v^//3) = (3h{x, t) . (22) 

Here h on the right side is defined by jSJ) with Z(x,t) as in item (A) and parameters 
v = ~, D = 7, and A = 0. 

The reader might ask, why we make such an extensive list. Firstly, the items 
provide a better understanding under which physical conditions the KPZ equation 
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is a valid approximation. Secondly, the list ensures that the natural but somewhat 
formal procedure in (A) and (B) properly captures our understanding based on 
discrete models. Thirdly, and perhaps most importantly, the various approximations 
do not only give mathematical sense to the KPZ equation but also provide a tool by 
which, at least in principle, some properties of the KPZ equation can be computed. 
In this respect, (A) can be used for a short time expansion. (C) has not been of help, 
yet. But a related discretization, where the directed polymer is placed on Z x R, 
has been used to compute the free energy [13]. The identity in (D) is the gateway 
to the replica method as used for spin glasses and other disordered systems with the 
special feature that the matrix element on the right side can be analysed through 
the Bethe ansatz. While the replicas forces one to work with divergent series, the 
method has turned out to be a powerful tool. Our exact solution is based on the 
approach (E). But at an intermediate stage we have to rely on a deep theory for the 
PASEP developed by Tracy and Widom in recent years [HJ l4"5l 135] . 

3 Exact solution of the KPZ equation with sharp 
wedge initial conditions 

The pdf of the random free energy log Z(x, t) with Z(x, t) as in Section 2 has been 
computed for every x, t. We use a Fredholm determinant formula by Tracy and 
Widom [35J, which in case of 0-1 step initial conditions provides the probability 
distribution for the position of the m-th particle at time t and thereby the probability 
distribution of h e of (l2~2~]) . While this is a convenient starting point, to actually 
establish the limit in ( 1221) still requires an intricate asymptotic analysis. 

To explain the structure of the exact solution we first recall the pdf for the Tracy- 
Widom distributed random variable £tw of dSJ). Like the pdf of h(x,t), it is defined 
in terms of a Fredholm determinant of a symmetric integral operator with kernel 
over R x R. We introduce the Airy kernel 



where Ai is the standard Airy function. The corresponding operator acting on 
L 2 (R,dx) is denoted by Kai- K A \ is a symmetric projection. We also introduce the 
projection onto the interval [s, oo), denoted by P s . P S K A \P S is trace class for every 
s > — oo. Then 



The Fredholm determinant on the right can be defined through the eigenvalues Xj(s), 




(23) 



Prob(£ TW <s) = det(l - P a K M P s ) . 



(24) 



j = 1,2, . . ., of P s K Ai P s as 



oo 



det(l - P s K Ai P s ) = J](l - A,-(s)) . 



(25) 
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In fact, it is more appropriate to think of P S K A \P S as a large matrix by evaluating 

as 

djp = Xs(xi)K Ai (xi,Xj)x s (xj) , i,j = 1, • • .,N, (26) 

with {xj,j = l,...,N} suitably chosen, e.g. equally spaced, base points, and Xs 
the indicator function of the interval [s, oo). Then 

det(l - P s K Ai P s ) * det(l - AW) , (27) 

for large N. The optimal choice of base points and error estimates are discussed in 

m- 

For the solution of the KPZ equation we obtain 

{\/2v)h(x, t) = -{x 2 /Avt) - ^( 7 *) 3 + 2 log a + 7t £t , (28) 

with j t = (a 4 ^) 1 / 3 , a = (2iy)~ 3 ^ 2 XD 1 ^ 2 . We note the self-similar flattening of the 
deterministic droplet shape proportional to t and a uniform shift, also proportional 
to t. The relative fluctuations are of order t 1//3 with a random amplitude £ t which 
is of order 1 but has a pdf changing in time. In (J7|) we anticipated the special case 
x = and t large. The pdf of £ t is given by 

/oo 
du-f t e^ s - u) exp[-e 7t(s - u) ] 
-oo 

x ( det(l - P u {B t - P Ai )P u ) - det(l - P u B t P u )) . (29) 
Here B t has the kernel 

/oo 
dw{\ - e~ 7tw )~ 1 Ai(a: + w)k\{y + w) . (30) 
-oo 

and P\i is an one-dimensional projection with integral kernel Ai(x)Ai(y). The first 
factor under the integral in (129]) is the pdf of the Gumbel distribution known from 
extreme statistics. (129]) looks like the convolution of two pdfs. However from nu- 
merical simulations we know that the second factor, while normalized, may take 
negative values. 

Equipped with the exact solution (|29[) . we can check whether its long time is 
given by the GUE Tracy- Widom distribution, to say whether £ t — > £ TW as t — > oo. 
Clearly, for t — > oo, the first factor of ( 1291) converges to 8(s — u) and the second 
factor to 

det(l - P u (K Ai - P M )P U ) - det(l - P u K Ai P u ) = — det(l - P u K Ai P u ) , (31) 

du 

the latter inequality following from identities proved in [47]. We conclude that the 
universal one-point pdf is established. For very short times has a Gaussian dis- 
tribution of width t 1 / 4 |39j. Hence our solution to the KPZ equation describes the 
crossover from a Gaussian at short times to the GUE Tracy- Widom at long times. 
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Plots of the pdfs covering mostly late times can be found in [49J. 

Note. The story behind the exact solution ( 1251) . ( 1291) is slightly convoluted. The 
formula was derived in the fall 2009 by Amir, Corwin, and Quastel [39] and inde- 
pendently by us (HI HH1 [50] . All papers use as their starting point a recent Fredholm 
determinant formula by Tracy and Widom for the PASEP with 0-1 step initial con- 
ditions valid for any q > 1/2, in particular for q = | + 2(3y/e [35]. Simultaneously 
the replica method was persued. Dotsenko and Klumov [5T1 E2] provide a detailed 
analysis of the eigenfunctions of the attractive 5-Bose gas on the line. With the 
hindsight from the exact solution, this leads to 

< exp [ - exp(-A + log Z(x, t) + (x 2 /Aut) + ±{ lt f - 2 log a)] ) 

= det(l - K x>t ) (32) 

for all real A and t > 0, expressing the generating function by a Fredholm determi- 
nant [531 G3] • The operator K\ t has the kernel 

K x , t (x,y) = (1 + e^ x -^)- l e^ x - X) K Ai (x,y). (33) 

Substituting A by 7 t A in (132|) and taking the limit t — > oo on both sides yields indeed 
the GUE Tracy- Widom distribution function for the rescaled height. Calabrese, Le 
Doussal, and Rosso [53] establish that if the average in ( 1321) is computed by using 
(129]) one indeed arrives at det(l — Kx,t)- (} 

Equipped with such an input the universality of the KPZ equation can be dis- 
cussed with more precision than before. One aspect concerns how well a microscopic 
system is approximated by the KPZ equation, which to be meaningful requires a 
tunable asymmetry. For example, the polynuclear growth model and also the liquid 
crystal of the experiment would have no such tunable parameter. But if available, 
then, roughly speaking, the description by the KPZ equation becomes valid in the 
limit of weak asymmetry and correspondingly long times. An example in case is the 
2D Ising model with Glauber dynamics at low temperatures with the H — interface 
oriented along the (1,1) direction. The asymmetry parameter is the external mag- 
netic field. For weak fields, which is needed anyhow in order to suppress nucleation 
in the bulk, the motion of the interface is well described by the KPZ equation. 

As a second aspect we note that the exact solution of the KPZ equation pro- 
vides us with qualitative information beyond the universal features. To illustrate 
we discuss the example of the approach to the Tracy- Widom distribution. In the 
experiment PQ, and also for the TASEP with 0-1 step initial conditions, one finds 
that the mean is the slowest mode, while higher cumulants decay rapidly to their 
Tracy- Widom value. Schematically, for large t we write 

(/i(0, t)) = Voo t + a.t 1 / 3 ^) + cot' 1 / 3 ) + Oir 1 / 3 ). (34) 

Here and a\ are model dependent parameters and we discuss the coefficient Cq. 
cq > for the experiment and TASEP. Relative to (£tw) the mean decays as i -1 / 3 
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with a positive amplitude. On the other hand one can expand ( |29l) in 1/t. The 
expansion is dominated by the Gumbel distribution, which produces a shift by t -1 / 3 
to the left. Thus the exponent comes out correctly, but the amplitude and the 
full first order finite time correction are model dependent. Since the KPZ equation 
corresponds to weak asymmetry, the PASEP with q close to 1/2 should have an 
approach also from the left. This is indeed confirmed by Monte Carlo simulations, 
which show a sign change of Cq at q w 0.7. 

4 Conclusions and outlook 

By our results we have added one prominent member to the KPZ universality class. 
Currently more universal pdfs have been computed for lattice models. For curved 
geometry and for flat initial conditions one knows not only the pdf for a single 
location but also all multi-point pdfs at the same long time, see [55J for the current 
status. In addition the case of stationary initial conditions has been studied. For 
the KPZ equation this would mean to choose an initial height such that dh(x, 0)/dx 
is white noise in x. One has to see whether the KPZ equation will catch up with 
these results in the future. 

From the point of view of the directed polymer, the continuum KPZ version is 
the only finite temperature model for which universal pdfs have been obtained so far. 
It would be of interest to understand how a comparable result can be accomplished 
for lattice models. 

Acknowledgements. We are grateful to Patrik Ferrari and Joachim Krug for most 
useful comments. 
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